Transcriptomics secondary analysis of severe human infection with SARS-CoV-2 identifies gene expression changes and predicts three transcriptional biomarkers in leukocytes.

SARS-CoV-2 is the causative agent of COVID-19, which has greatly affected human health since it first emerged. Defining the human factors and biomarkers that differentiate severe SARS-CoV-2 infection from mild infection has become of increasing interest to clinicians. To help address this need, we retrieved 269 public RNA-seq human transcriptome samples from GEO that had qualitative disease severity metadata. We then subjected these samples to a robust RNA-seq data processing workflow to calculate gene expression in PBMCs, whole blood, and leukocytes, as well as to predict transcriptional biomarkers in PBMCs and leukocytes. This process involved using Salmon for read mapping, edgeR to calculate significant differential expression levels, and gene ontology enrichment using Camera. We then performed a random forest machine learning analysis on the read counts data to identify genes that best classified samples based on the COVID-19 severity phenotype. This approach produced a ranked list of leukocyte genes based on their Gini values that includes TGFBI, TTYH2, and CD4, which are associated with both the immune response and inflammation. Our results show that these three genes can potentially classify samples with severe COVID-19 with accuracy of ∼88% and an area under the receiver operating characteristic curve of 92.6--indicating acceptable specificity and sensitivity. We expect that our findings can help contribute to the development of improved diagnostics that may aid in identifying severe COVID-19 cases, guide clinical treatment, and improve mortality rates.

Comparison of Intracellular Transcriptional Response of NHBE Cells to Infection with SARS-CoV-2 Washington and New York Strains.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported in Wuhan, China in December 2019 and caused a global pandemic resulting in millions of deaths and tens of millions of patients positive tests. While studies have shown a D614G mutation in the viral spike protein are more transmissible, the effects of this and other mutations on the host response, especially at the cellular level, are yet to be fully elucidated. In this experiment we infected normal human bronchial epithelial (NHBE) cells with the Washington (D614) strain or the New York (G614) strains of SARS-CoV-2. We generated RNA sequencing data at 6, 12, and 24 hours post-infection (hpi) to improve our understanding of how the intracellular host response differs between infections with these two strains. We analyzed these data with a bioinformatics pipeline that identifies differentially expressed genes (DEGs), enriched Gene Ontology (GO) terms and dysregulated signaling pathways. We detected over 2,000 DEGs, over 600 GO terms, and 29 affected pathways between the two infections. Many of these entities play a role in immune signaling and response. A comparison between strains and time points showed a higher similarity between matched time points than across different time points with the same strain in DEGs and affected pathways, but found more similarity between strains across different time points when looking at GO terms. A comparison of the affected pathways showed that the 24hpi samples of the New York strain were more similar to the 12hpi samples of the Washington strain, with a large number of pathways related to translation being inhibited in both strains. These results suggest that the various mutations contained in the genome of these two viral isolates may cause distinct effects on the host transcriptional response in infected host cells, especially relating to how quickly translation is dysregulated after infection. This comparison of the intracellular host response to infection with these two SARS-CoV-2 isolates suggest that some of the mechanisms associated with more severe disease from these viruses could include virus replication, metal ion usage, host translation shutoff, host transcript stability, and immune inhibition.

Analyzing the impact of a real-life outbreak simulator on pandemic mitigation: An epidemiological modeling study.

An app-based educational outbreak simulator, Operation Outbreak (OO), seeks to engage and educate participants to better respond to outbreaks. Here, we examine the utility of OO for understanding epidemiological dynamics. The OO app enables experience-based learning about outbreaks, spreading a virtual pathogen via Bluetooth among participating smartphones. Deployed at many colleges and in other settings, OO collects anonymized spatiotemporal data, including the time and duration of the contacts among participants of the simulation. We report the distribution, timing, duration, and connectedness of student social contacts at two university deployments and uncover cryptic transmission pathways through individuals' second-degree contacts. We then construct epidemiological models based on the OO-generated contact networks to predict the transmission pathways of hypothetical pathogens with varying reproductive numbers. Finally, we demonstrate that the granularity of OO data enables institutions to mitigate outbreaks by proactively and strategically testing and/or vaccinating individuals based on individual social interaction levels.

Impact of Local and Demographic Factors on Early COVID-19 Vaccine Hesitancy among Health Care Workers in New York City Public Hospitals.

Despite the development of several effective vaccines, SARS-CoV-2 continues to spread, causing serious illness among the unvaccinated. Healthcare professionals are trusted sources of information about vaccination, and therefore understanding the attitudes and beliefs of healthcare professionals regarding the vaccines is of utmost importance. We conducted a survey-based study to understand the factors affecting COVID-19 vaccine attitudes among health care professionals in NYC Health and Hospitals, at a time when the vaccine was new, and received 3759 responses. Machine learning and chi-square analyses were applied to determine the factors most predictive of vaccine hesitancy. Demographic factors, education, role at the hospital, perceptions of the pandemic itself, and location of work and residence were all found to significantly contribute to vaccine attitudes. Location of residence was examined for both borough and neighborhood, and was found to have a significant impact on vaccine receptivity. Interestingly, this borough-level data did not correspond to the number or severity of cases in the respective boroughs, indicating that local social or other influences likely have a substantial impact. Local and demographic factors should be strongly considered when preparing pro-vaccine messages or campaigns.

Preprocessing of Public RNA-Sequencing Datasets to Facilitate Downstream Analyses of Human Diseases

Publicly available RNA-sequencing (RNA-seq) data are a rich resource for elucidating the mechanisms of human disease;however, preprocessing these data requires considerable bioinformatic expertise and computational infrastructure. Analyzing multiple datasets with a consistent computational workflow increases the accuracy of downstream meta-analyses. This collection of datasets represents the human intracellular transcriptional response to disorders and diseases such as acute lymphoblastic leukemia (ALL), B-cell lymphomas, chronic obstructive pulmonary disease (COPD), colorectal cancer, lupus erythematosus;as well as infection with pathogens including Borrelia burgdorferi, hantavirus, influenza A virus, Middle East respiratory syndrome coronavirus (MERS-CoV), Streptococcus pneumoniae, respiratory syncytial virus (RSV), severe acute respiratory syndrome coronavirus (SARS-CoV), and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We calculated the statistically significant differentially expressed genes and Gene Ontology terms for all datasets. In addition, a subset of the datasets also includes results from splice variant analyses, intracellular signaling pathway enrichments as well as read mapping and quantification. All analyses were performed using well-established algorithms and are provided to facilitate future data mining activities, wet lab studies, and to accelerate collaboration and discovery.

Genome-wide bioinformatic analyses predict key host and viral factors in SARS-CoV-2 pathogenesis.

The novel betacoronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused a worldwide pandemic (COVID-19) after emerging in Wuhan, China. Here we analyzed public host and viral RNA sequencing data to better understand how SARS-CoV-2 interacts with human respiratory cells. We identified genes, isoforms and transposable element families that are specifically altered in SARS-CoV-2-infected respiratory cells. Well-known immunoregulatory genes including CSF2, IL32, IL-6 and SERPINA3 were differentially expressed, while immunoregulatory transposable element families were upregulated. We predicted conserved interactions between the SARS-CoV-2 genome and human RNA-binding proteins such as the heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1) and eukaryotic initiation factor 4 (eIF4b). We also identified a viral sequence variant with a statistically significant skew associated with age of infection, that may contribute to intracellular host-pathogen interactions. These findings can help identify host mechanisms that can be targeted by prophylactics and/or therapeutics to reduce the severity of COVID-19.